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Abstract 

We derive bounds on the Dirac Yukawa couplings of the neutrinos in seesaw models using the 
recent Large Hadron Collider (LHC) data on Higgs decays for the case where the Standard Model 
(SM) singlet heavy leptons needed for the seesaw mechanism have masses in the 100 GeV range. 
Such scenarios with large Yukawa couplings are natural in Inverse Seesaw (ISS) models since the 
small neutrino mass owes its origin to a small Majorana mass of a new set of singlet fermions. 
Large Yukawas with sub-TeV mass right-handed neutrinos are also possible for certain textures 
in Type-I seesaw models, so that the above bounds also apply to them. We find that the current 
Higgs data from the LHC can put bounds on both electron- and muon-type Yukawa couplings of 
order 10 -2 . 



I. INTRODUCTION 



Understanding the new physics behind small neutrino masses is a major current focus 
of particle theory research. The simplest scenarios for this are the various types of seesaw 
mechanisms, which have also formed the basis for understanding the observed mixing pat- 
tern, which is so different from the quark sector. Two of them, type-I [1] and the inverse 
seesaw [2] postulate the existence of heavy SM singlet fermions, generally denoted by the 
symbol N a ( a is a generation index) in addition to the SM fermions. The stability of 
these patterns can be guaranteed by extending the SM gauge group as we discuss below. 
These singlet fermions couple to familiar lepton doublets via Yukawa couplings of the form 
C-y = y a bL a HNb + h.c., which after electroweak symmetry breaking leads to the Dirac mass 
of the neutrino, which is expected to be in the GeV range. In the type-I seesaw, the tiny 
neutrino mass arises once we introduce an additional Majorana mass term M a bN a Nb for the 
TV's. Being SM singlets, these Majorana particles can take arbitrarily large mass values, 
which then guarantees that the neutrino mass given by the formula below is small: 



Note that since the scale of neutrino masses is known to be in the sub-eV range, the scale M 
is correlated with the magnitude of the Yukawa couplings y a b- For instance, if the N masses 
are in the 100 GeV range, it implies y ~ 10 -5 or so. However, for specific flavor pattern of 
both M a b and y a b, tiny neutrino masses can be realized without making y necessarily tiny. 
Examples of such theories are given in Ref. [3]. 

On the other hand, in the case of inverse seesaw models, in addition to the set of singlets 
N a , one adds another set of singlet fermions S a which form a Dirac mass M with the N a 
fields, i.e., M a bN a Sb- The model so far has conserved lepton number and leads to zero 
neutrino mass [4]. One then allows the S fields to have a Majorana mass matrix fig whose 
overall scale is in the keV range [2]. The neutrino mass formula in the case is given by: 



This also leads to tiny Majorana mass for neutrinos. The small scale of the us matrix can 
be explained in various ways [5] and is a window to further new physics. In this model since 
the small neutrino masses arise from the small values of //g matrix elements, the singlet 



M v = -yM- l y T vl 
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neutrinos (N, S) form a Dirac pair with masses in the 100 GeV range and yet have Yukawa 
couplings of order 0(1) without any fine-tuning. 

A key question for neutrino mass physics is how to test the seesaw mechanism [6]. The 
obvious first step would be to determine the couplings y a b and M a . The vast literature 
on neutrino mass models is devoted to precisely this question where additional theoretical 
assumptions such as symmetries are used for this purpose. The symmetries restrict the 
theory to a particular sub-space of the full parameter-space. If the seesaw scale is high, as 
is the case in most type-I models, this appears to be the only realistic possibility. We take 
a different approach in this paper. Since y a \, denotes the coupling of the iV-fermions to the 
Higgs boson, if M is in the 100 GeV range, experimental information on y can be obtained 
from the LHC data on the Higgs boson decays [7] . Our goal in this paper is to focus on this 
question. 

In this paper, we focus primarily on the inverse seesaw case and we will comment on the 
specific type-I scenarios where our results will apply with slight changes. 

This paper is organized as follows: in Section II, we discuss the minimal gauge models 
with inverse seesaw; in Section III we give an overview of the Higgs phenomenology of the 
seesaw models; in Section IV we shall illustrate how to constrain the parameter-space in 
seesaw models using the results of the Higgs searches in the final state llvv] in Section V 
we derive these constraints on the seesaw models by analyzing the rates of all the measured 
decay modes of the Higgs boson; the concluding remarks are in Section VI. In Appendix A, 
we list the relevant decay widths of the heavy neutrino. 

II. GAUGE MODELS WITH LARGE YUKAWA COUPLINGS AND SUB-TEV 
SINGLET FERMIONS 

There are two classes of models which can have 0(1) Dirac Yukawa couplings and ~ 100 
GeV mass for the singlet fermions responsible for the seesaw mechanism. One class, as 
noted above, is the inverse seesaw model, where this situation can be naturally realized, 
and a second class, which belongs to type-I seesaw-based neutrino mass models where with 
specific texture for the Dirac Yukawa matrices as well as for the singlet fermion mass matrix, 
one can have LHC-accessible parameters. Below we briefly discuss the two classes of models. 
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A. Inverse seesaw gauge model 



If we simply add two heavy SM singlet fermions iV a , S a to the SM, there are many gauge 
invariant terms in the potential and we will not get the simple inverse seesaw formula. The 
simplest extension of SM where the inverse seesaw formula arises naturally is SU(2)l x 
U(l)j 3R x U(1)b-l with gauge symmetry under which the singlet fermions N transform as 
(1, +1/2, —1/2) and S fields are singlets. The gauge symmetry is broken by a Higgs field 
x(l, +1/2, —1/2) acquiring vacuum expectation value (vev) and the SM-doublet H vev then 
breaks the gauge group down to electromagnetic U(l). The Yukawa Lagrangian for the 
lepton sector is given by: 



It is clear that once we substitute (H°) = t> wk = v/\/2 = 174 GeV and (x) = v BL , we get 
the inverse seesaw matrix [2] 



with the mass matrix for the neutrinos given by Eq. (2) with M = J'vbl (Note that y, f,^s 
are in general 3x3 matrices). At the renormalizable level, there are no other terms allowed 
in the Yukawa Lagrangian by the gauge symmetry and therefore the above inverse seesaw 
formula for neutrinos is stable under radiative corrections. 

Our goal is to find constraints on y as a function of M. There exists an extensive analysis 
on the constraints for light singlet fermions [8]. Although this analysis addresses only the 
constraints on Majorana neutrinos, some of their results apply to our case as well. However, 
the limits derived in Ref. [8] for M > 60 GeV or so are very weak. Furthermore constraints 
from neutrino-less double beta decay [9] derived on heavy sterile neutrinos do not apply to 
this case since in our model, the TV and S form a pseudo-Dirac pair and lepton number is 
almost exactly conserved. 

In order to use the LHC data to explore constraints on y and M in the 100 GeV range, 
we will assume that (i) vbl ^wk and (ii) the mass of Re(x°) is heavy compared to the SM 
Higgs boson so that neither the heavy gauge boson associated with (B — L)-symmetry nor 
the interactions of Re(x°) affect the Higgs boson decay modes we consider. 



C Y = y a bL a HN b + fabN a xS b + i^s ab S a S b + h.c. 



(3) 
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FIG. 1: The Higgs decay modes into 2d.2v mediated by the ISS couplings. 



It follows from the above Lagrangian that if one of the singlet fermions has mass in the 
100 GeV range, it will affect the Higgs branching ratios: for instance if < Mh, then this 
opens up a new mode for SM Higgs decay, i.e., h — > v a Nb, and the collider signal will arise 
from N — v mixing diagram in Fig. 1 where N — > uZ, £W. Folding W, Z decays, one will 
get final states with vvi a l^, where in the final state both charged leptons and anti-leptons 
will appear and the existing LHC data on these final states will provide constraints on y. 
Clearly, which charged lepton appears will depend on the flavor structure of y and /. For 
/ we will go to a basis so that it is diagonal i.e. the N are mass eigenstates with S field 
providing the chiral Dirac partner. 

B. Type-I seesaw case 

Turning to the type-I noted earlier, in generic models, the Dirac Yukawa couplings 

are very small for the seesaw scale in the TeV regime. However, for specific textures for y, 
it is possible to attain singlet fermion mass in the 100 GeV range with Dirac Yukawa y's 
of order 0(1) while still satisfying the neutrino oscillation data. In this case the singlet 
fermions could show up at the LHC. Two examples of this type of texture are [3]: 



with the constraint that jj^- + + = where M a are the singlet neutrino masses and 
the singlet fermion mass matrix is diagonal. Seesaw of this kind leads in the leading order 
to zero neutrino masses; small Sy a b can then generate the neutrino masses and mixings, and 
this will not affect the Yukawa couplings or the masses M a . 




(5) 
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A second example has y a b of the following type: 



Uab = 



Va $a e a 

Vb h £b 
\ Vc 8c e c J 



(7) 



with right handed neutrino mass matrix of the form 



/ Mi \ 



M = 



M 1 



(8) 



V M 2 J 



with e at b,c, 8a t b,c *C U a ,b,c- in the limit of e, 5 = the neutrino mass matrix is a null matrix 
regardless of the values of M a and y a ,b,c an d small neutrino masses arise from choosing e's 
and 5's small. Therefore again in this case, LHC data should put constraints on y a ,b,c vs M a . 
Since in this case, y's are arbitrary, they could be of democratic type which will mean that 
there are multiple flavor final states. 

Two important points regarding these models are worth emphasizing. In both models, 
as the small parameters (//g in the case of inverse seesaw, and 5 and e in the type-I model 
discussed above) go to zero, the neutrinos become massless. As a result, the neutral Higgs 
decays to two modes: h-^v + N^v + N. Unlike the case where both v and iV are Majorana 
fermions (as in generic type-I seesaw models), the Higgs decay rate to v + N final states 
becomes twice as large. Thus, the discussions of Higgs decay in both classes of models are 
similar except for the flavor richness in case (B) compared to case (A) . In case (B) model 
(II), for y a ~ y b ~ y c , all flavors couple with equal strength to a single fermion flavor iV a 
implying that final state signal becomes a combination of all three flavors. Similar situation 
also occurs for case (B), model I if a, /3 ~ 1. 

III. SEESAW HIGGS PHENOMENOLOGY 

As discussed above, a new Yukawa interaction in the leptonic sector can potentially alter 
the Higgs phenomenology in a dramatic way. In fact, in the models considered here, the 
Yukawa coupling of the neutrino can be sizable - in principle, much larger than the largest 
Yukawa coupling involved in the decay of a light SM Higgs boson, i.e., the one with the 
bottom quark which is of order 0(1CT 2 ). 
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Here we shall describe the observable consequences of the new Yukawa coupling. Focusing 
on the inverse seesaw case, it is apparent from the neutrino mass formula in Eq. (2) that 
neutrino mass fits will largely be dictated by the combination of the flavor structure of the 
Us matrix and the y-matrix. 

In quark-lepton unified theories, the natural expectation is that flavor mixings in y are 
"weak". So we can assume y to be a diagonal matrix. The next question is the relative 
magnitude of different flavor Yukawa couplings. Below we will consider the cases when (a) 
Vi> e Uva 3> y Vr , (b) y v ^> y Ve 3> y Ur , and (c) a democratic form where all the Yukawa 
couplings are equal. These different coupling structures imply distinctive flavor structure of 
the final states. From the interaction Lagrangian, Eq. (3), it follows that y a b mediates the 
decays illustrated in Fig. 1: 



where charged and neutral leptons arise in all the flavor combinations allowed by the form of 
the y matrix. The gauge bosons Z and W arise in the decay of N through the v — N mixing, 
that is ^=yvjm^ b . The gauge bosons in turn can decay into leptons or hadrons, such that 
the new Yukawa coupling contributes to the rates of both fully leptonic and semi-leptonic 
final states. 

The heavy particles Nj,, W and Z in the decay chain of Eq. (9) can be either on-shell or 
off-shell depending on the mass of iV 6 . In particular, when rrih > rriN b the Higgs decay width 
in Eq. (9) scales as y 2 while for < rriN b , it scales as y A . 

The fully leptonic decay mode of the Higgs llvv is among the most sensitive search 
channels of the Higgs and therefore is a very sensitive probe of the ^/-coupling, especially for 
rriN b not far above or below m^. This allows us to derive a bound on the (mjsr b ,y) parameter 
space as described in details in section IV. 

For definiteness in what follows, we shall focus on the case of the ISS model with only 
one light flavor of heavy neutrino which we choose to be the electron flavor, N e , and the 
corresponding Yukawa coupling y Ve to the Higgs doublet. We generically denote the mass of 
the heavy neutrino as mjv, suppressing the explicit flavor index e. 

In this scenario the decay chain of Eq. (9) becomes 



h -)> u a N b + c.c. -)■ v a v h Z + [u a £ c W + + c.c. 



(9) 




(10) 
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The decay of the Z boson to leptons produces pairs of charged leptons of all flavors, and 
hence, it contributes to the same flavor final states, namely 

h — > e~ e + u e u e and h — > /i~ ^v^v^ . 

The decay of the W boson produces a single charged lepton. In this case, the W is 
produced by a mixing effect of the heavy N e , and therefore, the W is always produced in 
association with an electron or a positron. Hence, the decay process mediated by the W 
only contributes to electron final states, namely 

h — > e~/i + z/ e z/ M + c.c. and h — > e~e + v e b , e . 

Altogether, due to the assumed coupling structure y Ve ^> y v ^> y Vr , the final state leptons 
of the Higgs decays mediated by the ISS are both of opposite flavor and of same flavor, with 
same flavor jj, getting its only contribution from processes with Z-boson intermediate state. 

For the other coupling patterns the resulting flavor structures are straightforward modi- 
fications of the one considered here. For the pattern (b) where y v ^> y Ue ^> y„ T , the main 
difference is that now it is the same flavor e final state that gets its only contribution from 
the Z-intermediate state. The opposite flavor final states again get contribution from both 
Z- and ^-mediated processes, though with the kinematics of e and /i exchanged. 

As in our analysis described in the next sections there is going to be little difference 
between e and /i, the bounds derived for the ISS with light N e also apply for the case of 
light N^. 

For democratic Yukawa couplings, we can assume a texture for the NS Dirac mass matrix 
so that the lightest heavy neutrino is a linear combination of all the three flavor states with 
equal probability, i.e., Ni = (N e + + N T )/y/3. In this case, if y is a unit matrix, then we 
expect all the possible flavor combinations to be equally populated as both the interactions 
of N a and the gauge bosons are flavor- universal. The analyses of the LHC experiments are 
mostly sensitive to e and fi flavors only, therefore for such a democratic case, only 4 out the 
9 possible final state flavor combinations are easily detectable x . Furthermore the effect of 

1 In principle, the leptonic decays of the r may be used to recover some sensitivity. However, they are only 
a fraction of the r decay modes, and there are detection inefficiencies. As such, the gain might not be 
huge. At any rate, the use of the r's requires careful experimental considerations that go beyond our 
competence; hence we prefer not to consider them. 
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the mixing of flavors in the lightest mass eigenstate N in general reduces the rate of e and 
\x leptons. Because of these effects the bounds should be recalculated carefully to account 
for the different flavor structure of the signal. 

However, neglecting the sub-dominant processes mediated by the Z, one can obtain an 
estimate for the bound in the democratic case. Indeed, the rates of the processes mediated 
by the W are just rescaled by mixing and multiplicity factors compared to those of the ISS 
model with only a light N e . As such, the bounds of the ISS case should be relaxed by a 
constant factor. In particular, from the fact that the rate of Eq. (10) scales as y% when 
rrih > rriN and as yf, e when rrih < m^, we find that the bound gets relaxed by a factor 
(3/2) x / 2 or (3/2) x / 4 in the two cases. 

In the following, we shall illustrate in details how to derive a bound on the Yukawa 
coupling in the ISS model. Using the same procedure, we also computed the bound for the 
type-I democratic case and we checked that this simple rescaling argument describes the 
actual bound to a very good accuracy. Hence, we give only the bounds for the ISS case, 
and those for the democratic type-I case can be obtained by the simple rescaling described 
above. 

Of course a more targeted analysis from the LHC experiment with flavor structure explic- 
itly taken into account could distinguish between these coupling structures. For instance, an 
imbalance between the e and n same flavor final states can distinguish case (a) and case (b), 
while case (c) would be favored if all the flavor combination appear to be equally populated. 

Irrespective of the flavor structure of the Yukawa interactions, the new Yukawa coupling 
mediates new decay modes of the Higgs, which contribute to its total width, and hence, in 
these seesaw models, the Higgs width is larger than in the SM: 

T/i = Tsm + r seesaw . 

The seesaw contribution to the the total width is necessarily model-dependent. For the case 
of the ISS with dominant coupling y Ve , the width 2 is given by 

where we have normalized the Yukawa couplings such that the Dirac mass terms of all 
fermions are m D = ^= for v ~ 246 GeV. In the case of the democratic type-I seesaw, as in 

2 This formula sums h —> vN and h —> vN 
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Eq. (7), we can take ydemo = y a = Vb = hence the Lagrangian essentially contains three 
interactions of equal strength for the Higgs decay, and the Higgs decay width is the same as 
in the ISS model once 2/demo = 2/iss/3 is taken. 

In Section V we shall discuss in details how to use the information on the measured rates 
of the several Higgs decay final states to put a bound on the seesaw coupling. 

IV. BOUNDS FROM THE SEARCH OF THE HIGGS BOSON IN THE ilvv FINAL 
STATE 

The ATLAS and CMS collaboration have both found evidence for a Higgs-like particle 
at around 125 GeV. The main evidence for the new particle comes from final states with 
resonant two photons or four leptons [10-12]. In addition to the 41 and 77 searches the 
LHC and TeVatron experiments searched for a SM Higgs boson in several other final states, 
including in the final state Ilvv [13-17]. 

The experiments presented cut-based as well as multivariate analyses to put a bound on 
the mass of the SM Higgs boson. The bounds from the multivariate analysis are generically 
(slightly) more stringent than the one obtained from the cut-based analysis. However the 
multivariate analysis cannot be easily reproduced with our means, and hence, we shall only 
use the cut-based analysis to derive our bound. 

Here we shall reinterpret the results of Ref. [13] to extract a bound on extra sources 
of Ilvv events. To do this, we shall repeat the cut-based analysis of Ref. [13] on event 
samples generated by Monte Carlo tools - matrix elements computed with Madgraph5 [18], 
showered and hadronized with PYTHIA6.4 [ ] and detector response parametrized by 
Delphes1.9 [20]. Hadrons have been clustered into jet with the anti-kT algorithm as 
implemented in FastJet2 [21, 22]. 

In Ref. [13], the CMS collaboration performed several analysis on the Ilvv sample col- 
lected in the year 2011 with the LHC running at 7 TeV center-of-mass energy. In particular 
they made a basic selection on the leptons, jets and missing energy of the events depending 
on the flavor of the final state leptons. The cuts for the opposite flavor (OF) and same flavor 
(SF) cases are reported in Table I. Then the analysis is specialized for specific values of the 
SM Higgs boson mass and further cuts are devised. These cuts are collected in the Tables 
II and III for the Higgs mass hypothesis of 120 GeV and 130 GeV, respectively. These 
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Baseline selection for all 


OF(e/i) 


SF(ja/j) 


SF(ee) 


Tin = 1, n e = 1 


rifj, = 2 


n e = 2 


\r] e \ < 2.5, |r/ M | < 2.4 


M < 2.4 


|r? e | < 2.5 


AR ee > 0.3 


AiV > 0.3 


Ai? ee > 0.4 


mET > 20 GeV 


mET>40 GeV 


ma > 12 GeV 


ma > 20 GeV 


prju > 45 GeV 



TABLE I: Baseline selection of the CMS cut-based analysis, as in Ref. [13] 



Analysis A 12 o (tailored for SM Higgs m h = 120 GeV) 



OF(efi) 



p T ,e 2 > 10 GeV 



SF(ee, /in) 



PT,e 2 > 15 GeV 



p T ,h > 20 GeV 



mu < 40 GeV 



Acj) < 115° 



m T ,iimET £ [80, 120] GeV 



TABLE II: Cuts added to those in Table I in the analysis devised in Ref. [13] for the hypothesis 
m h = 120 GeV. 



analyses, which we call t4 120 and v4 130 , are the most sensitive to new physics connected to 
the Higgs-like particle recently discovered, and hence, these are the only analyses we are 
going to repeat in order to extract our bound. 

To obtain a bound on the ISS model we impose that the events yield of the ISS after 
the cuts of both the analyses A 12 o and A 130 is not larger than the total number of events 
allowed by each of the analysis. This is done as follows: for each hypothesis for in the 
ISS model we compute the yield of events after the cuts of CMS. This depends on the mass 
mjv and coupling y Ue and it is given by 



niss(m N , y v 



L ■ o h 



e SM~ 



r(/i ->■ ww* -> iivv) 



r«M + r 



ISS 



T(h — > uN e + c.c. — » v e ljtkV 
Tsm + Tiss 

(12) 



11 



Analysis Aisq (tailored for SM Higgs = 130 GeV) 



p T ,e 2 > 10 GeV 



SF(ee, /j,fj,) 



p T ,i 2 > 15 GeV 



pr,i x > 25 GeV 



m u < 45 GeV 



< 90° 



m T ,UmET G [80, 125] GeV 



TABLE III: Cuts added to those in Table I in the analysis devised in Ref. 
m h = 130 GeV. 



] for the hypothesis 



where L = 4.6 fb _1 is the integrated luminosity used in the analysis, ah is the total Higgs 
production cross-section taken from Ref. [23], j and k are flavor indexes e, fi, and esM and Sjk 
are the efficiencies of the CMS selections for the decays mediated by the SM decay channel 
WW* and by decays of the ISS, respectively. 

As we just want to illustrate here how to obtain an upper bound on the Yukawa, we shall 
use only a few representative values of m^r, namely 60, 100, 140 and 200 GeV. For lighter itln 
the LHC searches tends to be rather ineffective. In fact, in the ISS signal both the charged 
leptons originate from N — > Ijl^v and therefore the invariant mass of the two leptons cannot 
exceed m^. The cuts m u > 12 (20) GeV for OF (SF) leptons, needed to reject leptons from 
QCD decays, would remove completely the ISS contribution for m N < 12 (20) GeV. Also, 
for light mjy the bounds from other experiments are more stringent [ ]. The power of the 
LHC is the sensitivity to around and above the mass of the Higgs, which improves 
significantly as compared to the reach of previous direct bounds. 

We remark that our computation of the total events yield neglects the possible interference 
between the ISS and the SM contributions, which in general is a small effect due to different 
flavor and Lorentz structure of the decays. 

The selection efficiencies esM and have been computed with showered events passed 
through Delphes1.9 and hadrons clustered into jets using the anti-kT algorithm with jet- 
cone radius parameter R = 0.5 , as done in the CMS analysis. The obtained efficiencies have 
been rescaled such as to reproduce the SM Higgs boson event yield in Table 3 of Ref. [13] 
for rrih = 120 and 130 GeV for the analysis v4 120 and A 130 respectively. The results of the 
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used chain of simulation codes are rather realistic, indeed the rescaling factor is almost flat 
w.r.t the cuts and the differences between the CMS and Delphes1.9 efficiencies are within 
20%. 

For the total width of the SM Higgs boson Tsm and the partial width T(h — » WW* — > 
llvv) we take the reference values of Ref. [23]. 

The width Tjss due to the decays mediated by y Ue and the partial widths T(h — > uN e + 
c.c. — > vjljlkv) have been computed with Madgraph5. These widths are the source of 
the dependence of n ISS on y Ue . In particular, they scale as yl when m h > m N , and as 
when rrih < ttl^, i.e., the heavy neutrino that mediates the decays is off-shell. Therefore, 
as N becomes heavier than the Higgs boson we expect the bound to quickly become less 
stringent. For completeness, we report our computation of the width Tigs and of the partial 
widths T(h — > vN e + c.c. — > v e ljikv) in Appendix A. 

To derive the bound we compute the maximal y Ue such that 

niss(m N ,yv e ) < ng 5 (A) 

where 77.95(^4) is the 95% CL limit on the number of events after the selection from the analysis 
under consideration. From Ref. [13] we extracted n 95 (Ai 2 o) = 55.9 and n 95 (Ai 30 ) = 78.6. 
We take as the final bound, which we denote as 2/^, 95, the most stringent one between the 
two bounds obtained for the analysis A120 and A130. It turns out that the bounds derived 
using the analysis A120 are stronger for all the cases we have considered. 

The obtained bound for a fixed rrih = 125 GeV as a function of mjv is shown in Figure 2. 
For the cases where m h > m^, we exclude y Ue ^ 0.01 while for m h < m N couplings y Ue <; 1 
are excluded. 

V. BOUNDS FROM THE OBSERVATION OF A HIGGS-LIKE PARTICLE AT 
THE LHC 

Evidence of a new particle has been observed in the 2011 and 2012 LHC data [10-12]. The 
region of phase space where the excesses are concentrated suggests that they are originated 
by a Higgs-like scalar particle with mass 125 GeV. Assuming that the new observed state is 
indeed the Higgs boson which is also involved in the seesaw, further bounds can be obtained 
from a global study of the properties of the new particle instead of just using the bound on 
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m h =\25 GeV 




80 100 120 140 160 180 200 
/TJ N [GeV] 



FIG. 2: Bound on y Ve for m/j = 125 GeV as a function of the mass of the heavy neutrino N e . 



the rate of new phenomena in the llvv channel. 

In what follows, we use the measured properties of the new particle and we shall illustrate 
how to use the ISS prediction to put a bound on the size of y Ue . 

The presence of additional decay modes h — > vN + Nv changes the properties of h in 
several respects. In fact the total width of h is increased w.r.t to the SM value. Deviations 
from the SM value of the total width are potentially observable in a line-shape analysis 
(when y Ue is large enough) or in a global analysis of Higgs decay data [24] . 

The change of the total width in turn affects all the branching fractions of h. In particular 
the rates of modes that do not get contributions from the new decay mode (such as 77, // 
and 41) are suppressed by a factor 

TlSS ^ . ^ ; 

J- SM + J- ISS 

where Tiss is given in Eq. (11). 

The Higgs decay mode h — > llvv instead gets a contribution from the decays in Eq. (10) 
and its rate is enhanced by a factor 

niss(m N ,y Ue ) 



(13) 



«ISS(^7V,0) 

where riiss is given in Eq. (12). 

Altogether the ISS model, compared to the SM, predicts a suppression by a factor 7 IS g in 
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the observed rates of all channels but llvv 3 which instead is enhanced by the new decays. 

The constraint that comes from the increase of the total width is far-reaching. In fact, it 
also applies to seesaw models where, for any reason, the mode llvv does not get enhanced. 
This might happen, for instance, when the N is very light, say itln < 20 GeV, such that the 
leptons from the N decay do not pass the selection cuts in Table I. Furthermore, the bound 
coming from the extra contribution to the width applies to other models as the case of ISS 
models where the coupling y Ur dominates the decay. In this case the searches into llvv 
are much less effective, still the effect on the total width provides a bound on the seesaw 
coupling. 

The constraint from the llvv channel are more specific to each model. Also they depend 
on the details of each of the WW analysis. In fact to make a correct use of the measured best- 
fit signal strengths of the WW channels one should consider each WW analysis separately 
and compute in details the efficiencies as we have done for the analysis of Ref. [13] discussed 
in Section IV . 

For these reasons we shall consider the bound from two sets of measurements. In both 
cases we take the ISS model as reference. In the first case we put a bound on y Ue coming from 
all the measured signal strengths excluding the WW channels. This is the safest possibility 
as we are making the least number of assumptions on the structure of the seesaw couplings. 
Additionally we compute the bound using all the available data on the Higgs-like particle at 
125 GeV. To deal with the several WW analyses we make the simplifying assumption that 
all the efficiencies for the processes mediated by the seesaw coupling are the same as those 
computed for the analysis of Section IV 4 . 

To put a bound we proceed as follows. For each searched decay mode of the Higgs the 
collider experiments give the best fit value of {j^pCp^j x BRd where p runs on the Higgs 
production modes and BR^s are the branching fraction for the various Higgs boson decay 
modes d. These measured best fits are expressed in units of the SM prediction and are 

3 Here we are disregarding the ISS enhancement of final states jjvv and Ivj j which, for a light Higgs boson, 
are less sensitive due to large backgrounds. 

4 While this is not completely rigorous we expect it to be a good approximation. In fact, the new CMS 
analysis [14] is very similar to Ref. [13]. Furthermore, one should note that the bound on the seesaw 
coupling is mostly sensitive to the ratio esfciAjfc in Eq. (12). As we find in section IV that the ISS 
efficiencies are quite similar to those of the SM, one can expect our simplifying assumption to be reliable 
for the ATLAS analysis as well. 
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CMS 77 2011+2012 


1.6 + 0.4 


CMS ZZ 2011+2012 


0.8 + 0.4 


CMS WW 2011+2012 


0.6 + 0.4 


ATLAS 77 2011+2012 


1.4 + 0.5 


ATLAS ZZ 2011+2012 


1.3 + 0.6 


ATLAS WW 2011 


0.6 + 0.6 


CMS bb AP 2011+2012 


0.1 + 0.6 


CMS rr 2011+2012 


-0.2 + 0.8 


ATLAS bb AP 2011 


0.5 + 2. 


ATLAS rr 2011 


0.2 + 1.8 


TeVatron bb AP 


2.1 + 0.7 


TeVatron WW 


0. ± 1. 


CMS WW AP 


-1.7 + 3.5 


CMS 77 Dijet 2011 


4.2 + 2. 


CMS 77 Dijet Tight 2012 


1.3 + 1.6 


CMS 77 Dijet Loose 2012 


-0.6 + 2. 



TABLE IV: Signal strength best-fits extracted from Refs. [10—12, 25—30]. AP stands for associated 
production. The best-fit for the AP and Dijet analyses are extracted from [29, 30] . We consistently 
take the best-fit signal strength for = 125 GeV throughout. 

referred to as best-fit signal strengths 

where by \i* d we mean the central value and by 8\i* d we mean the symmetrized la error on 
the best-fit of the channel d given by the experiments. 

To place a bound we confront the ISS prediction with the latest best-fit signal strengths 
given by the experiments at the LHC [10-12, 26-30] and the TeVatron [25]. For convenience 
of the reader the best-fit signal strengths used in our analysis are reported in Table IV. 

As in the ISS model the couplings of the Higgs to the fermions, gluons and to the gauge 
bosons are not modified the production cross-sections are the same as in the SM. This allows 
us to compute the changes in the rates from the changes in the branching fractions only. 
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FIG. 3: The x 2 a s function of y Ve in the ISS model for = 125 GeV and tun = 60, 100, 140 
and 200 GeV. The horizontal black dashed line corresponds to x 2 = Xmm + 1 • The colored solid 
lines are for the x 2 °f the entire dataset, and the dashed lines are for the x 2 without the WW 
channels. 

From the ISS predictions for the signal strengths fJ*d(yv e ) and the measured signal 
strengths we compute the x 2 f° r several choices of m N as a function of the ISS coupling 



where d runs on all the measured rates under consideration, and for the ISS \iii v9 is given 
by Eq. (13) and fid = liss for all the other channels. 

For each the \ 2 is minimized w.r.t. y Ve at a value Xmm- Considering all channels but 
those with h — > WW we get Xmm — 16, adding the WW channels we get xLia — 20. In all 
cases we find that the \ 2 is minimal for y Ue = 0, therefore we derive a 67% CL upper-bound 
on y v& . In Figure 3 we show the obtained x 2 — Xmm with an d without the WW channels. 
The inclusion of the WW channels improves the bound by a factor of about 2 in most cases. 



x 2 = E 



d 



MVuJ -/irf) 2 
(« 2 
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VI. DISCUSSIONS AND CONCLUSION 



In this paper we have derived bounds on the Dirac Yukawa couplings of the lepton doublet 
in inverse seesaw models for neutrino masses using LHC Higgs data. In generic versions of 
these models, such bounds are useful since one could understand small neutrino masses while 
keeping the Yukawa couplings to be of order one and the singlet fermion masses in the 100 
GeV range. We have focused on the cases where the electron or muon Yukawas are the 
dominant ones and also discussed the case with flavor- democratic Yukawa. Our discussion 
applies to the supersymmetric version of the model as well. It is perhaps worth pointing 
out that in SUSY ISS model, there are additional D-term contributions of order of a few 
GeV to the Higgs mass [31] thus relieving some MSSM parameter space. We find that for 
singlet fermion masses between 60-140 GeV, useful bounds can be derived on the Yukawa 
couplings from the recent LHC data on Higgs searches. 

It is also worth noting that in the case of inverse seesaw, there are limits on the mixing 
parameter || from leptonic unitarity [32]. The current bounds for the e-flavor is < 0.044 
and for the muon sector it is < 0.03. These bounds are weaker than what we obtain in this 
paper from LHC data. 
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Appendix A: Decay Widths of the Heavy Neutrino 

In this Appendix we collect the partial and total widths of the heavy neutrino N that 
have been used to compute the event yield in Section IV and V . Also we collect here the 
partial widths of the Higgs and the increment in the total width used in Section IV and V 
to compute our bounds. 
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m N 
(GeV) 


T{N e -> e~e + u e ) 


T(N e -> v e iTn + ) 
I^-(GeV) 


T(iV e -> e-/i+i/ M ) 
I^-(GeV) 


^•(GeV) 


60 


1.464 x 10" 4 


2.569 x 10~ 5 


2.32 x 10~ 4 


0.002716 


100 


0.03204 


0.001495 


0.03182 


0.3263 



TABLE V: The relevant partial widths and the total width of N e for m^r < m/j in the ISS model 
with dominant y Ue . 



(GeV) 


m N 
(GeV) 


Tiss 


Tiss(e + e ) 
[^•(GeV)] 


riss(/^> ) 

[^•(GeV)] 


riss(e =F / u ± ) 


125 


60 


2.9458 


0.1588 


0.0279 


0.2516 


100 


0.6446 


0.0633 


0.0030 


0.0629 



TABLE VI: The relevant widths h — > livv and the increment of the total decay width of the 
Higgs in the ISS model with dominant y Ue and < mh- The partial widths are calculated by 
multiplying the total width Tigs by the branching fractions of N computed from Table V. 



M h 
[GeV] 


M N 
[GeV] 


T(h -> u e W + e~) 
[V% • (GeV)] 


T(h -> v e Zv e ) 
[V% ■ (GeV)] 


Fiss(h) 
[V 4 N ■ (GeV)] 


125 


140 


1.658 x 10~ 4 


1.051 x 10~ 4 


5.42 x 10~ 4 


200 


1.119 x 10" 5 


6.655 x 10" 6 


3.57 x 10" 5 



TABLE VII: Total increment of the decay width of the Higgs in the ISS model with dominant y, 
for M N > M H . 



m h 
(GeV) 


m N 
(GeV) 


riss(e + e ) 
K ■ (GeV)] 


rissG«V ) 

K • (GeV)] 


riss(e =F ^ ± ) 
K • (GeV)] 


125 


140 


3.882 x 10" 5 


3.681 x 10" 6 


3.53 x 10" 5 


200 


2.657 x 10~ 6 


2.433 x 10" 7 


2.416 x 10~ 6 



TABLE VIII: The relevant decay widths h — > tlvv for > nih in the ISS model with dominant 
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